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Abstract — Multimedia streaming applications with stringent 
QoS requirements in 60GHz mmWave wireless personal area 
networks (WPANs) demand high rate and low latency data 
transfer as well as low service disruption. In this paper, we 
consider the problem of robust reflector placement in 60GHz 
WPANs. Reflectors relay traffic from transmitter devices to 
receiver devices facilitating i) the primary communication path 
for non-line-of-sight (NLOS) transceiver pairs, and ii) secondary 
(backup) communication path for line-of-sight (LOS) or NLOS 
transceiver pairs. We formulate the robust minimum reflector 
placement problem and the robust maximum utility reflector 
placement problem with the objective to minimize the number of 
reflectors deployed and maximize the network utility, respectively. 
Efficient algorithms are developed to solve both problems and 
have been shown to incur less service disruption in presence 
of moving subjects that may block the LOS paths in the 
environment. 



I. Introduction 

The millimeter wave (mmWave) band has attracted con- 
siderable commercial interests due to the advance in low- 
cost mmWave radio frequency integrated circuit design. The 
mmWave band provides 7GHz unlicensed spectrum resource 
at the center frequency of 6OGH2Q which would enable 
many high data rate applications like high definition streaming 
multimedia, high-speed kiosk data transfer and point-to-point 
terminal communication in data center (TJ-|[3), etc. 

In contrast to many existing RF technologies such as 
2.4GHz WiFi radios, mmWave radio has several unique 
physical characteristics. First, the propagation and attenuation 
loss are much more severe in the 60GHz band. It is shown 
that the free space path loss in 60GHz is more than 20dB 
larger than that in 5GHz. The oxygen absorption loss is 
as high as 5-30dB/km. Furthermore, the penetration loss is 
also much higher through typical building materials IT). As 
a result, line-of-sight (LOS) path is the predominant path 
for signal transmission, while signals along the second-order 
and higher-order reflection paths are highly attenuated and 
often negligible. Second, to combat such significant signal 
degradation, directional antenna technology is essential in 
mmWave devices. By using directional antenna on both the 
transmitter and the receiver sides, mmWave radios can obtain 
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significant gain in the received signal strength, while incurring 
negligible interference to/from other radios [4|-|6|. In this 
paper, we consider an mmWave wireless personal area network 
(WPANs) equipped with directional antenna on all devices. 

In addition to high bandwidth demands, multimedia appli- 
cations in 60GHz WPANs also have stringent requirements 
on service disruption (defined as the frequency and duration 
of time that network connectivity is not available), which can 
occur due to changes in channel conditions such as LOS link 
blockage by moving subjects in the space. This motivates us 
to employ reflectors for two purposes. First, reflectors can 
be used to relay traffic from transmitters to receivers that do 
not have LOS connectivity. Second, reflectors can provide a 
secondary (2-hop) path in case of blockage on the primary 
(direct) path. We consider active reflector for the ease of 
control of reflection direction and signal boost. 

In this paper, two reflector placement problems in 60GHz 
mmWave WPANs are investigated: robust minimum reflector 
placement (RMRP) that attempts to find the minimum number 
of reflectors and their best placements from a set of candidate 
locations with bandwidth and robustness constraints, and ro- 
bust maximum utility reflector placement (RMURP) that aims 
to maximize network utility given a fixed number of reflectors. 
Two vertex-disjoint (except for the endpoints) paths (one called 
the primary path, and the other the secondary path) are provi- 
sioned between each pair of transmitter and receiver. Conse- 
quently, seamless switching to the secondary path is facilitated 
in event of channel degradation or blockage on the primary 
path avoiding service disruption. Robustness is characterized 
by the D-norm uncertainty model, which models tolerance to 
concurrent (worst-case) failures of a subset of primary paths 
|7|. Using linear relaxation and the duality theorem of linear 
programming, RMRP and RMURP are transformed to the 
mixed integer linear programming (MILP) and mixed integer 
non-linear programming (MINLP) problems. Two algorithms, 
bisector search and Generalized Benders Decomposition, are 
proposed for RMURP and are shown to have near optimal and 
optimal performance. Extensive simulations demonstrate the 
fault tolerance of the proposed reflector placement algorithms 
in significantly reducing the probability of service disruption. 

The rest of this paper is organized as follows. The related 
literature work is analyzed in Section [II] In Section III we 
introduce the network model and the problem statements 
for robust reflector placement. The RMRP formulation is 
presented in Section [TV] while RMURP is formulated in a 
similar way in SectionjV] Furthermore, we also discuss two 
proposed algorithms for RMURP in Section [Vj Performance 
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evaluation is presented in Section VI Finally, we conclude 



this paper in Section VII 



II. Related Work 

Significant prior literature have been produced on differ- 
ent aspects of 60GHz radios, from CMOS circuit design to 
network protocol development. In this section, we summarize 
prior work on MAC design in mmWave WPANs. 

A spatial time-division multiple access (STDMA) scheme 
was proposed for a realistic multi-Gbps mmWave WPAN 
in p). With the help of a heuristic scheduling algorithm, 
it is able to achieve significant throughput enhancement as 
much as 100% compared to conventional TDMA schedules. 
In |9j, Cai et al. presented an efficient resource management 
framework based on the unique physical characteristics in a 
MC-DS-CDMA based mmWave networks. The authors also 
conducted extensive analysis of spatial multiplexing capacity 
in mmWave WPANs with directional antennae in fl0| , (TTJ. 
In J5;|, Madhow et al. conducted a probabilistic analysis of the 
interference in an mmWave network, as the result of uncoordi- 
nated transmission. It is concluded that an mmWave link can 
be abstracted as a "pseudo- wired link" with negligible interfer- 
ence when the beam width is 20 degree. Similar observations 
are made in Q, fl2) . Therefore, the primary interference at 
the transmitter or receiver devices is the predominant source of 
contention. In fl3) , 1 14 1, to address the deafness problem in- 
duced by directionality, Gong et al. proposed a new directional 
CSMA/CA protocol for IEEE 802.15.3c 60GHz WPANs. With 
virtual carrier sensing, the central coordinator can distribute 
the network allocation vector (NAV) information, to avoid 
collisions among the devices occupying the same channel. 
The author also extended the work to a multiple-user scenario 
in JT5|. A distributed scheduling protocol is proposed by 



coordinating mmWave mesh nodes in (6), and can achieve high 
resource utilization with time division multiplexing (TDM). 
However, none of the above work model or address reflector 
placement problems in mmWave WPANs with directional 
antenna. 

There are also some existing work on repeater selection and 
relay operation scheduling in mmWave WPANs. Repeater se- 
lection was investigated in 1 16 1, with the objective to maximize 
data rate for each transmitter and receiver pair by determining 
the best link allocation. In fTT) , Lan et al. explored time slot 
scheduling for relay operations in the scenario of directional 
antenna on mmWave devices and formulated the throughput 
maximization problem as an integer programming problem. 
However, both schemes do not consider robustness in presence 
of uncertain link blockage. 

To the best of our knowledge, we are the first to explore ro- 
bust reflector placement in 60GHz WPANs. Some preliminary 
results are presented in fl8). 



III. Network Model and Problem Statement 

Consider an mmWave network consisting of a set of L 
logical mmWave links (simplified as logical links), each link 
i € L is associated with a source device (transmitter) Sj, a 



destination device (receiver) dj, and a traffic demand bps. 
mmWave reflector devices (simplified as reflectors) equipped 
with steerable antennas can relay data between the transmitters 
and receivers. The reflectors can be placed at a set of K 
candidate locations. We further consider a set of O obstacles 
in the environment with known locations. 

A. Geometric Model for Link Connectivity 

In this section, we introduce a geometric model to char- 
acterize link connectivity in 60GHz mmWave WPANs. LOS 
transmissions are feasible between a transmitter and a receiver 
if and only if the direct path between them is unobstructed and 
their distance is less than a threshold d. 

To bound the end-to-end latency, at most two-hop paths 
(involving one reflector) are allowed between any transmitter- 
receiver pai^j Thus, the connectivity of any logical link is 
determined by the visibility regions of its end points defined 
as follows: 

Definition 1: (Visibility region) Given a 2-D plane of in- 
terest, any two points (a, b) are visible to each other if the line 
segment between them does not intersect with any obstacles 
and the distance of a and b is no more than d. The visibility 
region V(a) of a point a in the plane is the bounded shape 
consisting of all unobstructed points no more than distance d 
from a. 

Consider transmitter and receiver a and b, the (binary) 
connectivity of logical link (a,b) is thus characterized by: 



X(a,b) 



1, iff V(a) n V{b) ^ 
0, otherwise 



(1) 



If X(a, b) = 1, logical link (a, b) is feasible (directly or via a 
reflector in a and b's overlapping visibility region); otherwise, 
it is infeasible. For the rest of the paper, we only consider the 
set of feasible logical links given by: 



n = {i| A(« i ,d i ) = l,VieL}, 



(2) 



where Sj, di are the transmitter and receiver of the i-th logical 
link respectively. 

Fig. [T] illustrates the notion of visibility region. There are 
two mmWave end devices, DEVI and DEV2. The shadowed 
area corresponds to the overlapped visibility region between 
DEVI and DEV2, which is the candidate region for reflector 
placement for DEVI and DEV2. The existence of overlapped 
visibility region is the necessary and sufficient condition of 
the connectivity of logical links, in both LOS and non-LOS 
(NLOS) cases. 

The relationship between feasible logical links and candi- 
date reflector locations can be modeled as an undirected bi- 
partite graph G(Q,1C,E), where JC = 1, ...,K is the set 
of candidate reflector locations. An edge e = (I, k) exists 
between a logical link I € £1 and a candidate reflector location 
k E JC iff k is in V(si) D V(di), where s/ and di are the end 
points of /. Thus, the set of feasible logical links that can use 

2 The proposed solutions can be extended to cases where longer paths are 
allows. 
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Fig. 1: Visibility region and overlapped visibility region. 

a reflector placed at the candidate location k as relay is given 
by: 

n k = {i\ke v(si) n v(di),yi e f>},vfc e tc. (3) 

B. Needs for robustness 

Reflectors serve two purposes: i) providing the primary 
communication path for NLOS logical links; and ii) providing 
secondary (backup) communication path for LOS or NLOS 
logical links. Provisioning of secondary paths reduces service 
disruption when the primary path is obstructed. 

To see the impact of secondary paths, we conduct a sim- 
ple simulation study. Consider a home-network environment 
in Fig. |2ja), where there is a LOS logical link I and a 
dedicated reflector at a fixed location. In the robust setting, 
one reflector is used to provide a secondary communication 
path for /. Inside the room, there are M moving human 
subjects modeled as a circle with a radius of 0.3 meters. We 
adopt the random walk model [19], where in each step, a 
person moves 0.3 meters with the direction randomly chosen 
from the set {-90°, -45°, 0°, 45°, 90°}. Without reflectors, 
the communication between TX and RX is disrupted when 
a person blocks the direct LOS path. With reflectors, an 
outage occurs only when both the primary (LOS) path and 
the secondary path (via the reflector) are blocked. Fig. |2|b) 
and Fig. |2|c) show the percentage of blockage time and the 
mean, max and min duration with different number of human 
subjects, respectively. 

As shown in Fig. [2] (b)(c), when the number of moving 
subjects increases, the percentage of link blockage and block- 
age duration increase with and without the secondary path. 
However, the use of backup path reduces both the blocking 
probability and the duration of each outage. This translates to 
better quality of service (QoS) at the application layer. 

C. Problem Statement 

Reflector placement concerns with the selection of reflectors 
among a finite set of candidate locations to optimize for certain 
network utilities. We consider two variations of the problem. 

Definition 2: (Robust Minimum Reflector Placement 
(RMRP) problem) Given an mmWave network with a set 
of feasible logical links with fixed traffic demands, find the 



minimum number of reflectors and their locations among 
candidate location set JC that satisfy connectivity, bandwidth, 
and robustness constraints. 

Definition 3: (Robust Maximum Utility Reflector Placement 
(RMURP) problem) Given an mm Wave network with a set 
of feasible logical links with base traffic demands, find the 
placement of at most m reflectors among candidate location 
set JC such that the ratio of the achievable rates over the base 
rate is maximized subject to robustness constraints. 

We restrict relay of data to reflectors only. In the robust 
formulation, for each feasible logical link, two vertex-disjoint 
(except for the endpoints) communication paths are provi- 
sioned, one as primary path, and the other as secondary 
path. Both the primary and secondary paths between mm Wave 
transmitters and receivers cannot be more than 2-hops. If a 
reflector serves as relay for more than one logical links, time 
division medium access (TDMA) scheduling is adopted. The 
interference among concurrent transmissions is assumed to 
negligible due to the high directionality of mm Wave commu- 
nications, which is consistent with the measurement results 
reported in (6). Therefore, the main sources of contention 
arise from the half-duplex constraint and multiplexing at the 
reflector nodes. 

IV. Robust Minimum Reflector Placement (RMRP) 

In this section, we present the analytical form of the RMRP 
problem. The following notations are used: 

• Primary indicator: Xik — 1 if reflector k is selected 
by logical link i as its primary path reflector; otherwise, 

x lk = 0; 

• Secondary indicator: = 1 if reflector k is selected by 
logical link i as its secondary path reflector; otherwise, 

Vik = 0; 

• Selection indicator: z k = 1 if reflector k is selected by 
at least one logical link; otherwise, z k = 0. 

• NLOS indicator: r/j = 1 if logical link i does not have 
a LOS path; otherwise, rji — 0. 

If a logical link i has a LOS path, no reflector is needed for 
the primary path; otherwise, one reflector should be selected 
for the primary path, which should satisfy the following 
condition: 

K 

^x ik = T]i,Vi e tt. (4) 
fc=i 

On the other hand, at least one reflector is needed to 
facilitate the secondary path, that is, 

K 

5^ Vifc = i,v*Gn. (5) 

In addition, a reflector cannot be used for the primary path 
and the secondary path simultaneously. Therefore, we have: 

Xik + Vik < l>Vi € fljt,Vfe. (6) 

As mentioned before, for simplicity, we assume each re- 
flector has only one half-duplex transceiver. Therefore, the 
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Fig. 2: Link blockage with and without reflectors due to moving human subjects in a 10m x 10m room with a fixed mmWave 
TX/RX and a dedicated reflector, M human subjects moving randomly inside. 



transmission time of relaying a unit data of logical link i via 
reflector k is ^ 



Tik = 



(7) 



R Si ,k Rk,di 

where R Si .k, Rk,d t are the mmWave data bandwidth between 
the source and the reflector, and between the reflector and the 
destination, respectively. For the AWGN channel, they can be 
modeled as: 



and 



1 



Wl0 S |_- T- PnD{Sz . k ) 

(when D(k, d t ) > 6), 



Rk,di — 



1 



P t G t G r 



(when D(k, d t ) > 6), 



(when D(si,k) < 6) 



(8) 



(when D(k,di)<Q) 



(9) 

respectively, where W is the channel bandwidth in Hz, P t is 
the transmission power, P n is the noise floor level, G t , G r 
are transmitter and receiver antenna gains, 7 is the large-scale 
path loss index, D(a,b) is the distance from a to b, and is 
a constant threshold on communications radius. 



For a reflector k, the TDMA scheduling for the associated 
logical links should satisfy: 



(10) 



where fj is the traffic demand of logical link i, Tik is the unit 
data relay time of i via reflector k. The first term on the left 
side represents the percentage of reflector capacity occupied 
by all the logical links using this reflector as their primary 
path reflector. The second term represents the protection 
function for the set of logical links using this reflector as their 
secondary path reflector. A protection function <?&(.) measures 
the robustness of an mmWave WPAN. Its meaning will be 



To this end, the RMRP problem can be formally stated as: 



minimize 



(11) 



subject to Constraints |4} — \10\ 

variables x ik ,y ik ,z k 6 {0, 1}, Vi G fi, k = 1, . . . , K 

A. Reformulation under the D-norm uncertainty model 



Several uncertainty models have been proposed in literature, 
including General Polyhedron, D-norm, Ellipsoid, etc |7(. In 
this paper, we adopt the D-norm uncertainty model that is 
characterized by the following protection function: 



3fc(yk,r) 



max 

:S k <zn k ,\S k \=T, 



E 

ies k 



yik'i'ik' 



(12) 



Under the D-norm uncertainty model, among the set of logical 
links ilk that can relay through reflector k, at most Tk 
links will be blocked simultaneously on the primary path and 
consequently transmit on their secondary path via reflector k. 
The maximization gives the worst case traffic loads induced 
on the reflector. 

Two special cases are of particular interest. If Tk = \flk\, 
then gfc(yk, r) = X«eQ fc VikTink- This means all logical links 
in Slfc fail simultaneously. This corresponds to the maximum 
robustness. In this case, more reflectors may be needed. At the 
other extreme, if Tk = 0, no logical link is blocked. Fewer 
reflectors are in use. However, there is little fault tolerance 
in the resulting reflector placement. Denote p = Tk/\&k\ 
the robustness index p is a parameter to tradeoff between 
robustness and resource usage. . 



further explained in Section IV-A 



Under the above Z?-norm uncertainty model, (|T0j) can be 
rewritten as: 

7 riiXi k riTik + max } yikriT ik < Zk, Vfc, (13) 

.4^ s k :.s k cn k ,\s k \=r k ^ 

( fT3| is not directly tractable since it involves an inner- 
optimization in the protection function. The protection func- 
tion can be reformulated an integer linear programming prob- 



5 



lem as follows Q: 

max Y~] y ik riT ik s lk , 

{0<s lk <l }v ^ k 

S.t. ^ S ik < Tfc; 

s ife G {0,1}, Vz G O fe 



(14) 



Consider a linear relaxation of the above problem where 
Sih € [0, 1]. Due to the linearity of the constraints, the optimal 
solution occurs at the vertices of the feasibility region. Hence 
the optimal solution s* k must be either or 1, as there is no 
gap between the integer linear programming and the linear 
programming solutions. 

Taking the dual of the linear programming problem (14) , 
we have: 



min Ik^k 

{Pifc>0} Vie n i .,9fc>0 



s.t. q k 



i £ Q k 

Pik > yikr t n k , 



(15) 

Substituting <TT3J into ( fTT) , we can obtain the equivalent 
formulation of trie RMRP problem as a mixed integer linear 
programming (MILP) problem as follows: 



min 

x,y,z,p,q 



£ 



Zk 



s.t. 



^ riiX ik riT ik + q k T k + ^2 Pik < Zk, Vfc 
ien fc i£f2 fc (16) 

qk +p-ik > Viknnk, Vi g fi fc ,Vfc 

Constraints |4j([5](((6]l0([8j((2} 

variables x ik ,y ik ,z k G {0, l},p tk >0,q k >0 

B. Hardness of RMRP 

We prove in Appendix [A] that RMRP is NP-hard. 

The MILP problem in ( fTo") can be solved by any MILP 
solver. In our implementation, we adopt the MILP solver of 
the IBM optimization tool - CPLEX p0| . 

V. Robust Maximum Utility Reflector Placement 
(RMURP) 

In contrast to RMRP, which tries to minimize the number 
of reflectors, RMURP aims to maximize the total utility of an 
mm Wave WPAN given a fixed number of reflectors. 

Let be the base traffic demand on logical link i. We 
allow r,i be scaled up/down according to network constraints. 
That is, the actual data rate supported is wri, where a is a 
scaling parameter. This formulation is particularly relevant for 
transferring multimedia content that allows adaptive encoding. 
The objective of RMURP is henceforth to maximize the total 
utility U of the network, given by U = J2i ar i- 

The constraints of RMURP are similar to those of RMRP, 
except that the TDMA schedulability constraint ( flu) now 
becomes 



and additional cardinality constraints need to be included: 

k 

where m is the maximum number of reflectors to be used. 
To this end, the RMURP can be formalized as: 



maximize 

x,y,z,a 



(19) 



subject to Constraints Q(|5}((^0l[8jl|9jirr8j([T7| 

variables x ik ,y ik ,z k G {0, l},a > 0,V«,fe 

A. Reformulation under D-norm Uncertainty Model 

Again we c an apply the _D-norm uncertainty model (see 
Section HV-Al to the RMURP problem of ([19). This will 



Maximize 

x.y ,z,p,q,Q 



(17) 



transform (119 I into a mixed integer non-linear programming 
problem (MTNLP) as follows: 

E ar - 

i 

subject to 22 ViXikariTik + q k T k + ^ p ik < z k , Vfc 

isn fc ien k (20) 

q k +Pik > yi k ariT ik , Vi G Q. k ,Vk 
Constraints |4]l((5jl[6jl(7}l[8](l(9j([T8]l 

variables x ik ,y lk ,z k G {0,1}, p ik > 0,q k > 0,a > 

B. Hardness of RMURP 

We prove in Appendix [B] that RMURP is NP-hard. The in- 
clusion of variable a renders RMURP an MLNLP. Specialized 
algorithms need to be designed. 

Next, we propose two algorithms to solve the RMURP. 
The first algorithm is based on Bisection Search, which is 
shown to be fast but does not guarantee optimality. The second 
algorithm is based on the Generalized Benders ' Decomposition 
(GBD) technique pT[ , [22], which is proven to converge to 
optimal solutions but has higher computation complexity. 

C. Bisection Search 

Bisection Search is a heuristic method for finding roots of an 
equation. It iteratively bisects an interval and then selects the 
subinterval where a root must reside for the next iteration, until 
some termination condition is met. It is guaranteed to converge 
to a root of /(•) if and only if: / is a continuous function on 
the interval [A, B], and f(A) and f(B) have opposite signs. 

In ( |20] i, if a is given, RMURP becomes a MILP, which can 
be solved by CPLEX. The key is thus to determine the value 
of a. When a is large, RMURP is infeasible. When a = 0, 
RMURP is always feasible. More generally, if a,\ > ol-i, 
RMURP is infeasible when a — a\, then, RMURP is feasible 
when a = oi2- Treating feasibility and infeasibility as opposite 
signs, we apply the Bisection Search principle to decide the 
range of a iteratively until it is smaller than a threshold 
2 x TOL. Starting from an initial interval [A, B\, where 
a — A renders RMURP feasible and a = B renders RMURP 
infeasible, we substitute A or B with A ^ B depending on the 
feasibility of RMURP under a = A ~^ B . The "monotonicity" 
in the feasibility of RMURP with respective to a makes the 
Bisection Search converge fast but the optimality of the final 
results depends on TOL. 
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The Bisection Search based algorithm is summarized in 
Algorithm [T] 

Algorithm 1: Bisection Search 

Input : E2e data rate demands of every feasible mmWave logical link 
7"i,Vi; error tolerance TOL (which serves as the iteration 
termination condition); and the up-to-date known range for a: 
[A, B] 

Output: Maximum network utility Ut and reflector selection variables 
for every feasible mm Wave links x, y, z 

begin 

Set n = 1; 

while n < maxN do 

n , (A+B). 

L. < j , 

Solve the MILP problem RMURP(C). 
if RMURP(C) is feasible then 
A< — C; 

Obtain the solutions of RMURP (C): x( n ), y( n ), z( n ). 

else 

| B i — C; 
end 

n < — n + 1; 
if {B ~ A) < TOL then 

The best a found, ctbest < — ■A; 
Return U T , x( n ), y( n ), z< n ). 

end 

end 
end 



D. Generalized Benders' Decomposition (GBD) 

GBD is an iterative method for solving MINLP problems. 
The principle of the GBD algorithm is to decompose the 
original MINLP problem into a primal problem and a master 
problem, and then solve them iteratively. The primal problem 
corresponds to the original problem with fixed binary vari- 
ables. Solving the primal problem provides a lower bound, 
and Lagrange multipliers corresponding to the constraints. 
The master problem is derived through nonlinear duality 
theory using the Lagrange multipliers obtained from the primal 
problem. Solving the master problem provides a upper bound, 
and binary variables that can be used for the primal problem in 
the next iteration. It is proven to converge to the optimum 1 2 1 1 . 



Primal Problem: Let A := (x,y,z) represent the set of 
binary variables, A := (x,y, z) indicates the binary variables 
with specific values in {0, 1}. The primal problem 'P(A) of 
RMURP problem ( pO| i is obtained by fixing all the binary 
variables to A as follows: 

/(A) ^maximize ati 

p,q,ct — J 

i 

subject to ^2 rjiX ik ariTi k + q k r k + ^ p lk < z k , Vfe 

<?fc +Pik > VikariTik, Vi G fi fe ,Vfc 
variables P^r0,q>r0,a>0 

(21) 

Problem ( f2T| is a linear programming problem, which can 
be solved by any linear programming solver. Since the optimal 
solution of V(A) is also a feasible solution to ( |20] i, the optimal 
value /(A) provides a lower bound to the RMURP. It is also 
clear that, not all the choices of given binary variables can lead 



to a feasible primal problem. We need to treat it differently 
depending on whether the primal problem is feasible or not: 
• Feasible Primal: 

If the primal problem is feasible, let 

Tfc(A, p, q, a) = z k — ( ViZikOriTik + q k T k + ^ Pik), Vfc, 

g ifc (A, p,q, a) = q k + pik - y lk ar i T ik ,'ii G Q k ,\/k. 

(22) 

Then, we can compute the partial Lagrangian function for the 
primal problem as follows: 

L(A, P. q, X, v) = ^2 ari + XkTk + Uik &k , (23) 

i k k i 

where Xk,^ik > 0,Vi € Ofc,Vfc are the Lagrange multipliers. 
Thus, the Lagrange dual problem of V{A) can be stated as: 



min max L(A, p, q, a, X, v). 

\,v p,q,a 



(24) 



Since the problem is convex and has linearity constraints, the 
duality gap is 0. Thus, solving the Lagrange dual problem 
would give the optimal solution for V(A). 
• Infeasible Primal: 

If the primal problem is infeasible, we first define a set A 

as: 

A = {A\T k > 0,g ik > 0,Vi G fi fc ,Vfc,for some p,q,a}, (25) 
and consider the following feasibility-checking problem J 7 (A): 



minimize S 

P,1,a 



< 



subject to \2 ViXikanTik + q k Y k + Pik — z k < 5, Vfc 

yikarir ik - q k - pik < 5, Vi G Q fc , Vfc 
variables p^0,q^0,Q>0,<5>0 

(26) 

It is straightforward to see that, for any given A, V(A) is 
infeasible if and only if J-"(A) has a positive optimal value 

S* > 0. 

The Lagrangian function for .F(A) can be presented as: 

G(A, p, q, a, fi, <x) =^Mfc(^Z Vi^ikomnk + qk^k + ^ Pik 

k i£f2& i£S~2fc 

- z k ) + ^2 a ^{y t kar l T tk - q k - Pik), 

k i^Q.k 

V{^k,Cik) G 13 

(27) 

where fi/. , cr^ are Lagrange multipliers and 

13 = {{Hk,Gik)\v>k,o-i k > 0,X; fc (Mfe + E ie n fc CT «fe) = 
l,Vi€ fifc,VA:}. 

The Lagrangian dual of J 7 (A) becomes: 



max min G(A, p, q, a, /x, a). 



(28) 



Therefore, for any A G A, it can be characterized by the 
inequality constraint: 



> min G(A, p, q, a, fi, a). 

P,9 



(29) 
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Master Problem: The original problem in pO) can be written 

as: 



max 

A 



max f(A) 



= max 

AeA 



min max L(A, p, q, a, A, f) 

\,v p,q,a 



max /3 

s.t. ft < max L(A, p, q, a, A, i/), VA, v X 

p.q,a 



AG {0, l}n A, 



(30) 



where the second equality is due to J24j) because of the zero 
duality gap. Incorporating ( |2"9| ) into (pUjC we finally obtain the 
master problem A4(p, q, aTX v i /•*! trjas: 



max /9 

A 



S.t. 



/3 < max L(A, p, q, a, A, v), VA, ^ 
p,q,t* 

0>minG(A,p,q,a, ft( j),V(/i,(T)en (31) 
p,q,Q 

Constraints Q0@(7){8}l[9)£iJ) 

AG {0,1}, /3>0 

Note that, the master problem has two inner optimization 
problems as its constraints, which need to be considered for 
all A, v and p, a. This implies that the master problem has a 
very large number of constraints. In order to obtain a solvable 
mixed-integer linear programming problem, we employ the 
following relaxation for the master problem at iteration n as 
suggested by pT) : 

p < L(A n , p n , q n , a™, A", v n ) + V A^(.)(A - A n ), Vn G P fc 
> G(A n ,p n ,q n ,a n ,//> n ) + V aG(.)(A- A"),Vn G 

(32) 

where "P fc and J 7 * 1 are the sets of feasible and infeasible primal 
problems solved up to iteration k, respectively. 

The relaxed problem provides the upper bound of the 
original master problem and also provide the value of binary 
variables for the primal problem in the next iteration. 
The GBD algorithm is summarized in Algorithm [2] 
It has been be proved in [21] that, the solutions to discrete 
variables A 1 ,...,A fe do not repeat. Therefore, due to the 
finiteness of the discrete variable set, GBD algorithm con- 
verges within a finite number of iterations. When it converges, 
the lower bound is equal to the upper bound. Thus, optimality 
is achieved. 

VI. Performance Evaluation 

In this section, we evaluate the performance of the reflector 
placement solutions using simulations. In the simulations, an 
mm Wave home network is deployed in a lOmx 10m room, 
where N mmWave devices and O obstacles are uniformly 
placed. The reflectors can be placed at any grid point in a 
grid separated by distance do. The transmission radii of all 
mmWave end devices and reflectors are set to 6 meters. The 
(base) traffic demand of each logical link i is chosen as | of 
the AWGN Shannon channel capacity of the slowest LOS path. 
The PHY parameters are given in Table [I] In all experiments, 
d Q = 2m, O = 10. 



Algorithm 2: Generalized Benders' Decomposition 

Input : Base logical data rate for each feasible logical link i 
Output: Maximum network utility adapt factor a and reflector selection 
variables for every feasible logical link A = (x, y, z) and p, q 

begin 

set n = 1 and choose A £ {0, 1}, 

LB < oo, UB° < — oo, P° < — 0, T° < — 0. 

while LB"" 1 < UB 71 - 1 do 

if the primal problem is feasible then 

Solve the primal problem r P(K n ) to obtain optimal 
solution p n , q n , a n 
and Lagrangian multipliers A™ , v n ; 
V n < — V™- 1 U {n}, T n i — J 7 "- 1 ; 
LB n < — max(LB n ~ 1 ,f(A n )); 
if LB n == f(A") then 

| (A*,p*,q*,5*) i — (A n ,p n ,q n ,5"); 
end 

else if the primal problem is infeasible then 

Solve the feasibility-check problem -F(A) to obtain the 
optimal solution p n , q" , a n and Lagrangian multipliers 

At n ,cr n ; 

■pn < pn ) jrn < -pn-1 y | n |. 

end 

Solve the master problem M (p" , q n , a n , A™ , u n , fi n , a n ) 
and obtain the optimal solution A n+1 and f} n \ 
UB n i — /3 n ,n < — n + 1; 

end 

return A*, p*, q*, a*. 
end 



TABLE I: PHY Parameters 



PHY parameters 


Values 


Channel 


AWGN with gain 1 


Path Loss 


free space, exponent 2 


Transmission power 


20mW (13dBm) 


Noise floor 


-lOOdBm 



A. RMRP Performance 

In this section, we examine the performance of RMRP under 
different configurations by varying the number of mmWave 
logical links (N), the number of moving subjects (M) and the 
robustness index (p). 

Fig. |3ja) shows the number of reflectors when reflector 
candidate locations and obstacles are fixed, while the number 
of mmWave logical links varies. In this set of experiments, all 
logical links are feasible. More reflectors are needed as the 
number of logical links increases. However, the relationship 
is not always linear due to the absence of LOS paths between 
TX/RX pairs and the multiplexing of reflectors. 

Fig. |3jb) show the percentage of blockage time per link 
when human subjects move randomly in the room. The mo- 



bility setup is similar to that in Section III-B Clearly, As 
the number of human subjects increases, the percentage of 
blockage time increases as well. However, the robust scheme 
leads to 50% less blockage. 

Next, we evaluate the impact of robustness index p. In this 
setup, 1 human subject moves randomly and there are 5 logical 
links. Figure [3|c) shows the number of reflectors used and the 
percentage of blockage time. As expected, as p increases, more 
reflectors are used and the blockage time reduces. 
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# of reflector used 
Percentage of blockage time 
O : 




# of mmWave logical links: N 

(a) Effect of the number of logical links 
with p = 1 

Fig. 3: Performance of RMRP in an mm Wave home network deployed in a lOmxlOm room. 



# of moving subjects M 

(b) Effect of the number of moving 
subjects with p = 1, N = 5 



Robustness index p 

(c) Effect of the robustness index with 

M = 1,N = 5 



B. RMURP Performance 

We now evaluate the performance of the Bisection Search 
and GBD algorithms on RMURP. The error tolerance of 
Bisection Search is set to TOL = 1.0. 

Fig. [4] shows the utility achieved using both methods by 
varying the number of logical links (Fig. |4j a )X the total 
number of reflectors (Fig. |4jb)) and the robustness index 
(Fig.^c)). In all cases, GBD achieves higher utility compared 
to the Bisection Search method. Reducing the threshold TOL 
improves the performance of the Bisection Search but comes 
at a higher computation cost. 

Fig. |5ja) demonstrates the convergence of the GBD algo- 
rithm. As shown in Fig. |5Ja), over time, the upper bound 
(solutions to the master problem) is non-increasing; and the 
lower bound (solutions to the primary problem) is non- 
decreasing. The algorithm converges to the optimal solution 
after 50 iterations when the upper bound equals to the lower 
bound. 

Fig. |5jb)(c) shows the percentage of link blockage time per 
link under RMURP when human subjects move randomly in 
the room by varying the number of moving subjects and the 
robustness index, respectively. In both cases, GBD achieves 
lower percentage of blockage. This implies that the reflector 
selection in GBD has more spatial diversity. 

VII. Conclusion 

In this paper, we formulated robust reflector placement 
problems in mmWave WPANs, namely, the robust minimum 
reflector placement problem (RMRP) and robust maximum 
utility reflector placement (RMURP) for better connectivity 
and robustness against link blockage. Under the D-norm 
uncertainty model, RMRP and RMURP were casted as MILP 
and MLNLP problems. Efficient algorithms were devised and 
evaluated using simulations. 

For future work, we will incorporate more complex inter- 
ference and antenna models in the robustness formulation. 
Another interesting agenda is to explore the use of passive 
reflectors in mmWave WPANs. 
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Appendix 

A. Proof of the NP-hardness of RMRP 

Definition 4: We call the special case of RMRP problem, 
where the robustness index p = 0, the MRP problem. 

Lemma A.l: MRP is NP-hard. 

Proof: Since p = => T k = in MRP, the protection 
function of ( fT3| l is null, which means only primary paths are 
considered in ( flO) . In MRP, represents the percentage 

of reflector fc's capacity occupied by logical link i's primary 
path, if i choose to route its primary path via k. If we treat a 
bin as a reflector's capacity, and treat the volume of an item 
as the percentage of reflector capacity occupied by an logical 
link, the Bin-Packing problem can be reduced to MRP in a 
polynomial time. Since Bin-Packing is NP-hard |23|, MRP is 
also a NP-hard problem. ■ 

As MRP is just a special case of RMRP, RMRP is hence- 
forth harder than MRP, so RMRP is also NP-Hard. 

B. Proof of the NP-hardness of RMURP 

Definition 5: We call the special case of RMURP problem, 
where the mmWave network adopts robustness index of p = 
and only has one candidate reflector location, as MURP 
problem. 

Lemma A.l: MURP is NP-hard. 

Proof: Since p = 0^r fe = 0in MURP, the protection 
function term is null. Also, consider only one candidate 
reflector location k in the network, the scheduling constraint 
of MURP in ( fl7| i can be rewritten as: ^\ rjiXikCtriTi < z^. 

Let ari represent the value of item i, rjiariTi represent the 
weight of item i. the 0-1 knapsack problem can be reduced to 
MURP in a polynomial time. Since the Knapsack problem is 
NP-hard J23), MURP is NP-hard. ■ 
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Fig. 4: Performance of two RMURP solutions in an mmWave home network deployed in a lOmx 10m room, Bisection Search 
algorithm with 1.0 tolerance and GBD algorithm. 
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Fig. 5: Performance of RMURP in an mmWave home network deployed in a lOmxlOm room. 



As MURP is a special case of RMURP, RMURP is hence- 
forth harder than MURP. This proves the lemma. 
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